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ABSTRACT 

Wc combine recent estimates of dust extinction at z ~ 4 — 7 with UV luminosity function (LF) 
determinations to derive star formation rate (SFR) functions at z ~4, 5, 6 and 7. SFR functions 
provide a more physical description of galaxy build-up at high redshift and allow for direct comparisons 
to SFRs at lower redshifts determined by a variety of techniques. Our SFR functions are derived from 
well-established z ~ 4 — 7 UV LFs, UV-continuum slope trends with redshift and luminosity, and 
IRX-/3 relations. They are well-described by Schechter relations. Wc extend the comparison baseline 
for SFR functions to z ~ 2 by considering recent determinations of the Ha and mid-IR luminosity 
functions. The low-end slopes of the SFR functions are flatter than for the UV LFs, Aa ~ -1-0.13, 
and show no clear evolution with cosmic time (z ~ — 7). In addition, we find that the characteristic 
value SFR* from the Schechter fit to SFR function exhibits consistent, and substantial, linear growth 
as a function of redshift from 5 MQyi:~^ at z ~ 8, 650 Myr after the Big Bang, to ^ 100 M0yr~^ at 
z ~ 2, ~ 2.5 Gyr later. Recent results at z ^ 10, close to the onset of galaxy formation, are consistent 
with this trend. The uniformity of this evolution is even greater than seen in the UV LF over the 
redshift range z ~ 2 — 8, providing validation for our dust corrections. These results provide strong 
evidence that galaxies build up uniformly over the first 3 Gyr of cosmic time. 

Subject headings: Galaxies: high-redshift — Galaxies: evolution 



1. INTRODUCTION 



In recent years great progress has been made 
in determining the UV luminos i ty fun c tion (LF) a t 
early times (e.g., iSteide l et al.l. (l999l: iWvder et al 



2005t lArnouts et a l.. 20051 : iReddv fc Steidell. 120091 : 
Bouwens et all . I2007i . i2011lf lOesch et al.l I20ia r2012bl ) 



Recent results on the UV LF indicate a bright end that 
builds up substantially with time, with the character- 
istic luminosit y increasing from the highest redshifts to 
z ^ 3 — 4 (e.g. iBouwens et al.l . [2Q07D . Furthermore, the 
faint-end slope of t he UV LF is found to be very steep at 
redshifts z > 3 (e.g. IBouwens et all . [20Trb[ ) . This bright- 
ening of and steep faint-end slope a are consistent 
with a general picture where galaxies build up hierarchi- 
cally and where lower luminosity galaxies play a major 
role in the reionization of the universe. 

Despite the general usefulness of the UV LF for prob- 
ing early galaxy formation, one significant limitation of 
the UV LF in this regard is the sensitivity of rest-frame 
UV light to dust extinction. Inferring this attenuation 
directly at z > 3 is challenging since only the most bolo- 
metrically luminous systems at these redshifts can be 
detected in far-IR. Also, for these very high redshifts, 
tracers of star formation such as X-ray, radio, 24^m and 
Ha are either too faint to observe or redshifted out of 
the observable wavelength range of the most sensitive 
telescopes. 
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Fortunately, we can take advantage of a relation 
betw een the UV -c ontinu um slope /3, with fx cx 
A'' (jMeurer et all . Il999f) . and the likely dust ex- 
tinction to convert t he observed U V-lum inosities 
into SFRs (see e.g. IBouwens et"an . I2009D . It 
has only now been possible to accurately establish 
the distribution of the UV-continuum slope distri- 
bution ov er a wide ran g e in lumin o sity and red- 
shift (e.g. IBouwens et a .1. I2009L 120121: IWilkins et al 



20111: iFink elstein et all. 120111: iCastellano et al.l . 12012 



DunloD et a l.. 2012). This is the result of the installation 



of the WFC3 /IR camera on the Hubble Space Telescope 
and the recent acquisition of dee p wide-area dat a over 
the HUDF09+CA1 N 1DELS fields (IBouwens et all . [20T1 : 
iGrogin et all l25TlHKoekemoer et al.Ll2011| ). 

In this paper, we utilize these new UV-continuum slope 
determinations to derive SFR functions at z ~ 4, 5, 6 
and 7. SFR functions can better connect the growth 
of galaxies at the highest redshifts to galaxies found at 
later epochs, as well as give insight into the manner in 
which star formation lights up dark matter halos. Previ- 
ously, such functi ons were presented a t z ~ — 1 based 
on UV-hIR dat a (iMartin et all . [20051 : iBell et all . [20071: 
iBothwell etHl . \2QlM- We begin this paper by describ- 
ing how we use /3 to produce dust-corrected UV LFs and 
SFR functions (§2). In §3 we compare our SFR functions 
with the literature over a range of redshifts (z — 8). 
For ease of comparison with previous studies we adopt 
Ho = 7 0kms-^Mpc-\ = 0.3 and 17a = 0.7. We 
adopt a iSalpeteil jl955f) IMF with limits 0.1-125 
throughout this paper. 

2. SFR FUNCTIONS 

This section describes how we compute SFR functions 
from published UV LFs and UV-continuum slopes /?. We 
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Figure 1. Top: The relation bet ween tlie UV-continuu m slope /3 
and UV luminosity at 2 ~ 4 from iBouwens et al.l (201% with the 
binned points and linear fit < /3 > = — 0.11 (Myv ab + 19-5) — 
2.00. The inset panel shows the luminosity de-trended distribution 
of UV-continuum slopes around the mean relation. Middle: The 
average dust extinction from the /3-distr ibution in the upper panel 
as a function of luminosity, assuming a IMeurer et al.l I I1999I ) dust 
correction, as described in i|2.1l The errorbar in the bottom left 
corner shows both the random and systematic uncertainties in the 
relation. Bottom: The black and red points show the z A UV 
LF fr omlBouwens et al.l 1)20071 ') before and after correction for dust 
(see ^2.111 . The errorbar in the top left corner shows the fiducial 
error in the average dust correction. The dust-corrected UV LF 
has a fiatter faint-end slope a and brighter M^y-. 

begin by deriving the SFR functions in a fully stepwise 
fashion, applying our luminosity-dependent dust correc- 
tions to binned UV LFs (Section 2.1). We then derive 
analytical formulae for these SFR functions based on the 
UV LFs we use as inputs (Section 2.2). 

2.1. Dust-corrected luminosity functions 

We correct UV LFs for the effects of dust attenuation 
using the well-known correlation of extinction with the 
UV-continu um slope P. W e take the IRX-/3 relation es- 
tabhshcd bv IMeurer et all (|I999), 

Ai6oo = 4.43 + 1.99/3. (1) 

IMeurer et al.l (jl999l ) estimated the relation based on 
starburst galaxies in the local universe. Similar re- 
lations have been found at z ~ by other groups 



Table 1 

Stepwise determinations of the SFR fu nctio n at 2 ~ 4, 2 ' 
2 ~ 6 and z ~7 (see ^TJl 

logioSFR(MQyr-i) 0sFR (Mpc-'^dex-i) ' 



—0.66 


0.05920±0. 02855 


-0.44 


0.06703±0.00838 


-0.21 


0.02537±0.00326 


0.02 


0.02534±0.00268 


0.25 


0.01430±0.00144 


0.48 


0.01153±0.00117 


0.72 


0.00601±0.00025 


0.95 


0.00354±0.00017 


1.19 


0.00221±0.00012 


1.44 


0.00139±0.00008 


1.68 


0.00052±0.00006 


1.92 


0.00023±0.00004 


2.16 


0.00002±0.00002 


2 ~ 5 


-0.33 


0.01766±0.00858 


-0.11 


0.01161±0.00294 


0.12 


0.01076±0.00121 


0.36 


0.00420±0.00046 


0.61 


0.00362±0.00040 


0.86 


0.00224±0.00014 


1.11 


0.00121±0.00008 


1.37 


0.00060±0.00006 


1.63 


0.00023±0.00002 


1.89 


0.00006±0.00002 


2 ~ 6 


-0.04 


0.01197±0.00262 


0.41 


0.00426±0.00089 


0.77 


0.00173±0.00037 


1.01 


0.00110±0.00024 


1.26 


0.00026±0.00008 


1.51 


0.00014±0.00004 


1.77 


0.00002±0.00002 


2.03 


0.00002±0.00002 


2 ~ 7 


-0.07 


0.01543±0.00473 


0.15 


0.00761±0.00215 


0.38 


0.00513±0.00149 


0.61 


0.00224±0.00075 


0.84 


0.00106±0.00044 


1.08 


0.00031±0.00019 


1.32 


0.00033±0.00018 



(|Burgarella et al.l . [20051: lOverzier et all l20lT[) Though 
there is some evidence that the Me urer et al.l ()1999() re- 
lation does not work well for all sources, e. g. very young 
galax ies (< 100 Myr) and ULIRGs fe.g. iReddv et all . 
|2006( ). this relation has been found to be accurate in 
the mea n out to z ^ 2, despite consi derable scat- 
ter (e.g. iDaddi et al.l. l2007t iReddv et al.l [200 6. 20lfl 
Reddv et al.l . 120121: iMagdis et all . I2010al: iMagdis et all . 
There is even evidence in the new Herschel obser- 



vations that the IMeurer et al.l ()1999t ) relation is reason- 
ably accurate in estimat ing the dust extinct ion for 2: ~ 2 
Lyman-Break Galaxies (jReddv et al.l I2012D . We there- 
fore quite reasonably utilize this relation in interpreting 
higher redshift samples. 

Recently there have been a number o f studies exam- 
ining /? in high-redsh if t sam pl es fe.g. IBouwens et aL . 
200a IWilkins et all. [MTI: IBouwens et all . ~|20T2 



Finkelstein et al.l . 120111: lOunlop et a l.l l2012D. Perhaps 
the most definitive of these studies is IBouwens et al.l 
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Figure 2. Left: Both the analytical and stepwise SFR functions derived in this study from dust-corre cted UV LFs. Th e ste pwise SFR 
functions (individual points) were derived using the UV LF results f rom [Bouwens et al.l II2007I . I2011bl) . as described in ^2.11 The solid 
lines are the SFR functions derived in Schechter form as described in i|2.2l with parameters listed in Table |2l The lines are not fits to the 
points. The excellent agreement between the two approaches provides a useful cross-check. We have indicated the SFRs (>~ 150 MQyr~^: 
equivalent to a bolometric luminosity > 10^^ Lq) where we expect our SFR functions to be more uncertain due to incompleteness in the 
U V selections and possible unreliability of the I RK-fS relation (e.g. IReddv et al.1, [20061). The best e stimates at high SFRs may come from 
searches in the mid-IR/far-IR (black point from lDaddi et all . 120091 . see also lMarchesini et al.L [20101) . The SFR function therefore may fall 
off m ore slowly than we i nfer (dashed line). Right: A comparison of the SF R functions with similar functions derived from the bolometric 
LF of IReddv et al.l I I2008I : grey region), the IR LF of lMaenelli et al.l 112011k black open squares) and the Ha LF from lSobral et al.l I I2012I : 
black dashed line) at z ~ 2. The trend in the SFR, function, derived from our dust-corrected UV LFs at z ~ 4 — 7, clearly extends to z ~ 2. 
The smooth evolution in the SFR function provides some corroboration for the dust corrections we apply. 

to the UV luminosities of individual bins of the LF. We 
shift each point in the LF toward brighter magnitudes 
and correct for the fact that the luminosity-dependent 
dust correction increases the width of the bins. 

The bottom panel of Figure [1] shows the effect of our 
luminosity-depen dent dust correct i on on the stepwise UV 
LF at z 4 from iBouwens et al.l (j2007D . The dust cor- 
rection shifts the LF to higher luminosities, particularly 
at the bright end, causing M^y to brighten, and also 
flattens the faint-end slope. There is also a small shift to 
lower volume densities due to the rcnormalisation of the 
magnitude bins. 

Now that we have dust-corrected UV fluxes we can use 
well-established relations to compute the SFR as a func- 
tion of luminosity, giving us our d esired SFR functi ons. 
Wc use the following relation from iKennicuttl (fTOOl . 



(|20T1 . who using the CANDELS-fHUDF09 datasets 
find that /3 correlates with both redshift and luminosity, 
with higher redshift a.nd lower lumiiiosity g alaxies being 
bluer. The results bv IBouwens et al.l (|2012[ ) are in excel- 
lent agreem e nt with other result s in th e literature (e.g . 
Ouchi et all. I2004t iLabbe et all . I2007t lOverzier etall . 



20081: IWilkins et aT.I . l201lir 



For our dust corrections we assume a linear relation 
between the UV -continuum slope P an d luminosity, such 
as that given in IBouwens et al.l ()2012D . This is shown in 
the top panel of Figure [T] for z ~ 4, 



<I3> = 



df3 



dM 
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UV,AB 



+ 19.5) + 13 



A/uv = -19-5, 



(2) 



where 



dp 



dM„ 



and /3 



A/u 



-19.5 



arc from Table 5 of 



IBouwens et all (|2012[ ). Note that for z - 7 we use a 
fit with a fixed slope ,ff = —0.13 obtained from our 

^ rfMuv 

z ~ 4 — 6 samples, given the large uncertainties in this 
slope at z 7 and the lack of evidence for evolution in the 
/3-luminosity relation over the redshift range z ~ 2 — 6. 
The distribution of UV-continuum slope /3 shows sub- 
stantial scatter about relation 2, that can be approxi- 
mated by a normal distribution (Figure [TJ but remark 
there is a fatter tail toward redder /?). 

The steps for computing an average dust correction 
are as follows. We use Eq. [l]to calculate the UV absorp- 
tion ^1600 for each source in our adopted /3-distribution. 
Then to obtain the extinction correction at a given Myv 
we integrate over the /3-distribution (normal distribu- 
tion with mean < (3 >, Eq. [21 and cr^=0.34), setting 
^1600 ~ when Aiqqo < 0. The middle panel of Fig- 
ure [T] shows the resulting < Amuy > as a function of 
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(3) 



Since this relation gives the time-averaged SFR over a 
~ 100 Myr time window, it will underestimate the SFR 
(typically by 5, 2x) in ga l axies substantially younger 
than this (e. g . IVerma et al.l . I2007t IBouwens et all . 120121: 
IReddv et ahl . I2012D . However. Eq. |3] should work on 
average for the extended SF histories expected in LBGs. 

The left panel of Figure[2]shows the stepwise SFR func- 
tions at z ~ 4, 5, 6 and 7, based o n ste pwise UV LFs 
derived from IBouwens et al.l ()2007l ) and IBouwens et al.l 
(l2011bD . For convenience, we include our stepwise SFR 
functions in Table [T] 

2.2. Analytical SFR functions 
In this section, we use an analytical Schcchter-like ap- 



luminosity a.t z ^ A. We then apply this dust correction proximation to represent SFR functions at z 4 — 7. 
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Figure 3. Derived Sciiechter parameters agpR (top), (^sfr C™*'^" 
die) and SFR* [bottom) as a function of redshift. The red points 
are the derived dust-corrected SFR function parameters, while 
the black poi nts and crosses show e quivale nt p arameters for the 
UV L Ps from IBouwens etaT . ( 2007, ' 201 Ibl ') and IReddv fc Steidell 
II2009I ), respectively, before dust correction (black points are shifted 
+0.15 on the x-axis for clarity). The errorbars for our parame- 
ters include both the random and systematic uncertainties in 
and /3a/jjv=— 19-5 ■ compare to SFR parameters obtained from 
IBothwell cFaD (12011'; red diamonds), IBell e t al.l (20071; red plus 



signs), iReddy ct al. (2008; red triang les), [MaKnclli ct al. (2011; 
red asterisks) and Sobral et al.l l|2012l : red crosses), as described 
in jjS] Our dust corrections result in flatt er faint-end slopes a 
(Aa ~ +0.13) relative to the UV LP ( ^3^211 . Our dust correction 
also doubles the rate at which logio SFR* grows with cosmic time 
over the redshift range 2 ~ 4 — 7 ( ^3. 11 1. Clearly, dust correc- 
tions can have a significant impact on the apparent evolution of 
SF galaxies with time. 

We assume that the IRX-/3 relation for individual 
galaxies is described by Auv = Cq + Ci f3 and the distri- 
bution of galaxies at a certain A/yv is given by a Gaus- 
sian with fj,/3 =< P > (Eq. [5]) and width ap, which gives 

< ^M„v >= Co + 0.2 In 10 Cf(jl+Ci<l3>. (4) 
This expression is only valid in the limit that the dis- 



Table 2 

Schechcr parameters determined for the present SFR functions 



(^) 




^°SlO Mpc-3 


"SFR 


3.8 


1.54±0.10 


1.07+0.17 


-1.60+0.07 


5.0 


1.36+0.12 


0.76+0.23 


-1.50+0.12 


5.9 


1.07+0.17 


1.08+0.39 


-1.57+0.22 


6.8 


1.00+0.30 


0.64+0.56 


-1.96+0.35 



INotes. these bchechter param eters are obtained by dust 
correcting the UV LF using the lMeurer et al.l I) 19991 ) IRX-/3 
relationship. We adopt the linear relation between the 
UV-continuum slope /3 and UV luminosity (Eq. [2]l found by 
IBouwens eFaLl l l20l3) . Sec ^Ol 

tribution of UV-continuum slopes /3 does not extend 
to /5 ;$ —2.3 since such blue /3's formally give nega- 
tive dtist_correcMons_(a clearly non-physica l result) using 
the iMeurer et all (|1999f ) relation. For the iMeurer et all 
(| 19991 ) relation Eq. H simplifies to < Am^v >= 4.43 + 
1.82a|-hl.99</3>. 

To compute SFR function s we start with t he UV LF, 
described in Schechter form (jSchechterl . 119760 : 



dL^cj)* — 



Substituting SFR for L (and SFR* 
El and H yields 



L 
D 

for L' 



), using Eq. 



(/>(SFR) dSFR = 



'-^l dM 



X exp 



/ SFR 

lySFR* 

SFR 



dSFR 



SFR* / SFR* 



(6) 



where we have made the simplifying assumption that 
the cut-off in the Schechter function is exponential and 
not some slightly shallower high-end cut-off (the modi- 
fied functional form is consistent with the observations). 
This gives the conversions 



asFR 



)rr -f^l 



Ci 



df) 
dM 



UV ,uncorr 



1 



'-^l dM 



(7) 



(8) 



We calculate SFR* using Eq. [21 E] and HI i.e. we use 
logioSFR* = -0.4(M*v„r- < Ami,. >) - 7.25. 
Though Eq. [4] is a reasonable approximation for < 
^Muv >i our estimate of SFR* is slightly more accu- 
rate when we estimate < ^m*^ > as described in i i2.1l 
by setting Aigoo = when ^leoo < 0. Therefore wc will 
use this more accurate estimate of SFR* quoted in Table 

m 

The left panel of Figure O compares the analytical 
SFR functions (solid lines) with the corrected stepwise 
UV LFs described in ij2.1l They are in excellent agree- 
ment, providing a useful check on the analytic description 
used here. We would expect our derived SFR functions 
to be more uncertain at high bolometric luminosities 
(Lboi > lO^^L©: indicated on Figur e [H) where dust cor- 
rections are likely less reliable (e.g. iRcddv et al.l . 120061) 
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Table 3 

Schechter parameters for SFR functions in the literature (see also Figure [Sj 



(z) log 



SFR* 
1" AJq yr-l 



1„„ <PSFR 



"SFR 



probe 



reference 









0.91±0.01 


-3 


80±1 


09 


-1.51±0.08 


MIR+UV 


Bothwell et al. 2011" 





2 


0.88±0.27 


—3 


OlitO 


27 


— 1.45(fixed) 


MIR+UV 


Bell et al. 2007 





4 


0.97±0.15 


-2 


97±0 


06 


-1.45(fixed) 


MIR+UV 


Bell et al. 2007 





6 


1.06±0.12 


-2 


68±0 


10 


-1.45(fixed) 


MIR+UV 


Bell et al. 2007 





8 


1.19±0.15 


-2 


77±0 


15 


-1.45(fixed) 


MIR+UV 


Bell et al. 2007 





8 


1.10±0.07 


-2 


47±0 


08 


-1.56±0.12 


Ha 


Sobral et al. 2012 


1 


5 


1.41±0.06 


-2 


61±0 


09 


-1.62±0.25 


Ha 


Sobral et al. 2012 


2 


2 


l.TlitO.lO 


-2 


73±0 


12 


-1.57±0.23 


Ha 


Sobral et al. 2012 


1 


5 


2.28±0.15 


-3 


44±0 


23 


-1.60(fixed) 


MIR 


Magnelli et al. 2011 


2 





2.27±0.15 


-3 


41±0 


26 


-1.60(fixed) 


MIR 


Magnelli et al. 2011 


2 


3 


2.35±0.15 


-3 


49±0 


24 


-1.71±0.08 


MIR+UV 


Reddy et al. 2008 


3 


8 


1.54± 0.10 


-2 


97±0 


07 


-1.60±0.07 


UV+/3 


This work 


5 





1.36± 0.12 


-3 


12±0 


13 


-1.50±0.12 


UV+/3 


This work 


5 


9 


1.07± 0.17 


-2 


97±0 


16 


-1.57±0.22 


UV+/3 


This work 


6 


8 


l.OOit 0.30 


-3 


20±0 


38 


-1.96±0.35 


UV+/3 


This work 



theiKe 



attl 



Note s. These Schechter parameters are derived from the published Ha, bolometric, and UV LFs in the literature usinj 
||T998|) relations. 

° It is unclear i f there is a typographical error in the Schechter parameters provided by IBothwell et al.l II2011I ) . A simple integration of 
I 2011 ) SFR function does not yield the quoted SFR density. However we quote the numbers as they are presented in 
\ 2^T% . converted to a Salpetcr IMF with limits 0.1-125 M©. 



Bothwell et al. 



Bothwell et al 



and our Lyman-break selections may be more incomplete 
to dusty star-forininR ga laxies (e.g. iDaddi et al.l 120091 : 
iMichalowski et all 120101 ). While this should not affect 
the turnover in the SFR function, the SFR function may 
fall off less steeply than the exponential form adopted 
here. 

The analytical Schechter parameters arc presented 
in Table [2] The SFR*, 0gpj^ and asFR are calcu- 
lated assuming the z ~ 4 — 7 UV LF par ameters from 
iBouwens etaU (|2007| ) and iBouwens etall ()2011bl ) 

dMl. 



/3. 



Muv = -19.5 



and 



the 

0.34 from IBouwens et al.l 

(|2012D and the IMeurer et all (|1999( ) relation (Eq. H]). 
For the z ~ 7 LF parameters, we modify our procedure 
somewhat due to the fact that Eq. |3]is not especially ac- 
curate for the very blue /3's found at z 7 (sec discussion 
following Eq. H]) . Specifically, we keep the SFR* fixed to 
the value described above and then fit the low-end slope 
and the normalisation of the Schechter function to the 
stepwise UV LF. 

We have verified through Monte-Carlo simulations that 
we can successfully recover (AM* < 0.13 mag, Aa < 
0.03) the observed UV LFs and /3-luminosity relation us- 
ing the de rived SFR fu nctions and a /3-SFR relation (e.g. 
ICastcllan o et al.l . l2012[ ). so our approach can be inverted. 
We describe these simulations in detail in Appendix A. 

We remark that i f P shows a w e aker dependence 
on luminosity than IBouwens et al. (120 121) find (e.g. 

it would re- 
higher 



and 



iCastellano et al.l . 120121 : iDunlop et aL 
suit in a steeper asFR (Aa 0.1) anu mgnci v^gpR 
(A0gpp{ 0.05 dcx). Uncertainties in both the incom- 
pleteness and contamination corrections used for the 
UV LFs, of course, a.lso affect our SFR f unctions (e.g. 
IBouwens etldl . 120071: iGrazian et al] . I2011D . 

3. RESULTS 

One of the main results of this paper is our determi- 
nation of the approximate Schechter parameters for the 
SFR function. These parameters allow for an intuitive 
and quantitative look at how this function evolves with 



Table 4 

Indicative SFRs over the redshift range z < 



' 2 to 2 ~ 8 



<z) log 



A/q yr- 

(fixed t/)*)'^ 



(fixed n(> SFR))' 



probe 



2.0 


2.04±0.08 


2.01±0.30'-' 


MIR 


Magnelli et al. 2011 


2.2 


1.89±0.02 


1.80±0.30'= 


Hq 


Sobral et al. 2012 


2.3 


2.08±0.03 


2.05±0.30'^ 


MIR+UV 


Reddy et al. 2008 


3.8 


1.58±0.10 


1.52±0.30'^ 


UV+/3 


This work 


B.O 


1.33±0.11 


1.27±0.30'^ 


vv+p 


This work 


5.9 


1.09±0.16 


1.05±0.30'^ 


UV+^ 


This work 


6.8 


0.93±0.29 


0.86±0.30^ 


vv+p 


This work 


8.0 


0.61±0.04 


0.49±0.30^ 


UV 


Oesch et al. 2012a 



Notes. These Schechter parameters are derived from the 
published Ha, bolometric luminosit y and UV lumino sity 
functions in the literature using the lKennicutH 1)19981 1 relations. 
" Fixed logio 4'sf'R. ~ ~3-05. See also Figure |4] 

Fixed number densit y n(> SFR) = 2 • lO-^Mpc-^ 
llPapovich et al.l . [2011) . 

Indicative err ors from IPapovic h et al. l 1120111). I f we use a 
IChabrici' lf2Q03) IMF as in lPapovidT'et al .1 1|20111 ') (instead of a 
[Safpeter 1955: IMF with limits 0.1-125 M©) the SFRs given here 
would be 0.2 dex lower. 

cosmic time, as well as allowing for a comparison with 
lower redshift determinations. We discuss both the high- 
SFR and low-SFR end of the SFR function. 

3.1. High-SFR end: linear build-up of log SFR* from 
z - 8 to z - 2 

The high-SFR end of the SFR function is interesting 
since it provides us with a direct measure of the rate at 
which luminous galaxies are building up at early times. 
Wc compare our parameters with various studies that 
combine UV and IR ob servations to obtain complete SFR 
functi ons a t z < 2, i.e. IBothwell et al.l (|201lD . iBell et al.l 
(|2007[) and lReddv et all (1200 81 and also with Ha a nd IR 
LFs at z - 2 from lSobral et al . (2012) and Magnell i et al.l 
(|2011f ) respectively. We expect the IR LFs to probe 
the SFRs in the dominant population of dusty galaxies 
(ULIRG-|-LIRGs ) sX z ^ 2. We use the relations from 
iKennicufrB (|1998[ ) to convert the Ha and bolometric lumi- 
nosities to SFRs. The SFR functions at z ~ 2 are shown 
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Figure 4. The turnover SFR* at fixed 0gpj^ as a function of redshift (left) and cosmic time [right). The SFR* from lMagnelli et al.l II2011I ). 
ISobral et al.l II2012I ). [Reddy et al.l 1 )20081) and lOesch et al] I l2012al ) are indicated with a red asterisk, cross, tria ngle and square, respectively. 
The dashed line gives a linear fit to our derived logio SFR* at z ~ 4 - 7, with logio SFR* oc -0.22 (1 + z) ( ^JJl . The fit at z ~ 4 - 7 is 
in broad agreement with the literature over the entire range z ~ 8 to z ~ 2. The extrapolated fit also agrees with the current estimate 
for SFR* at z ~ 10 (jBouwens et al.L I2011al : lOesch et al.LI2012bl ). but the uncertainties in this estimate are still quite large. This provides 
strong evidence that galaxies build up consistently with time, from z ~ 8 to z ~ 2. 

in the right panel of Figure [21 these SFR functions are 
in reasonable agreement with each other, except at the 
high-SFR end where the Ha LF is low. This may result 
from an incomplete sampling (or inadequate dust eor- 
rections) of dusty galaxies by the Ha study. The z ^ 2 
SFR functions are consistent with the evolution observed 
between 2 ~ 4 — 7. 

The bottom panel of Figure |3] (see also Table [3]) shows 
SFR* for our results at z - 4-7 with the SFR* 
from pu blished studies mentioned a,bove. The black 
symb ols (|Bouwens et al.l . 120071 l2011bt iReddv fc Steidel 
|2009( ) represent the equivalent SFR* one would derive 
without applying a dust correction. There is a signifi- 
cant difference in the evolution of the high-SFR end of 
the SFR function taking dust attenuation into account: 
SFR* (without a dust correction) peaks at ^ 2OM0yr~^ 
for z ^ 3 — 4, while SFR* derived from dust-corrected 
UV or MIR LFs continues to increase from ^ 10 Moyr~^ 
at z ~ 7 to ~ lOOMeyr-i at z 2. 

In Figure |4] we show the evolution of SFR* at fixed 
<^SFR (with logio 0gFR = —3.05: see also Table |4]). This 
is interesting since it allows us to examine the evolution 
in the high-end of the SFR function without introduc- 
ing additional "noise" from the SFR*-0gpj^ degeneracy. 
We extend our eoin parison of SFR* to z ^ 8 usi i ig the 
Oesch et al.l (l2012al) UV LF (see also lTrenti eran . l20TTI : 



Bradlev et al.l . l2012D . The effect of the dust correction on 
SFR* decreases strongly with redshift for z > 4 and we 
expect the contribution at z 8 to be nearly negligible. 
We fit a linear slope to our own estimates of logio SFR* 
at z - 4 - 7 and find logio SFR* = 2.43 - 0.22 (1 + z) 
(Figure HI . The SFR functions without a dust correction 
follow logio SFR* oc —0.11 (1 -I- z). Comparing the two, 
we see that our dust corrections double the rate at which 
the high-SFR end of the SFR function grows with cos- 



mic time (see also the discussion in §6.2 of lBouwens et al.l 



[2003) . 

Extrapolating the best fit SFR* -redshift relation to 
higher and lower redshift we find good agreement with 
z 2 and z ~ 8 determinatio ns of SFR* ( Figure [j]). 
The estimate for z ~ 10 from lOesch et all (|2012bf ) is 
also consistent with the trend, but the uncertainty in 
SFR* at that redshift is still quite large. The evolution 
of SFR* with redshift at fixed 0* suggests that galax- 
ies build up in a consistent way during the first 3 Gyr 
of cosmic time. Moreover, the excellent agreement be- 
tween different probes of star formation across cosmic 
time, provides a valuable cross-check on the validity of 
t he dust correctio n s we u se. 

iPapovich et al.l (|2011[ ) also consider an evolving 
SFR vs. redshift relation based on an abundanee- 
matched selec t ion. Not surprisingly, given that 
iPapoyich et al.l (120111) use essentia lly the same UV LFs 
(from lBouwens et al.l . l2007l l2011bf ) and s imilar dust cor- 
rectio ns at z > 4 as we use here (fr om iBouwens et al.l 
120091 instead of IBouwens et al.l l2012f ) , they find similar 
SFRs over the range z ~ 4 — 8 (at constant number den- 
sity n{> SFR) = 2 X 10-''Mpe"^ see Table[4|). However, 
the present analysis suggests much more clearly that the 
build-up in the SFRs of galaxies extends not simply from 
z 8 to z 4, but continues all the way to z 2. 

3.2. Low-SFR end: flatter slopes and large uncertainties 

The low-end slope of the SFR function is interest- 
ing, since it provides us with insight into the physics 
and feedback effects in the lowest luminosity galaxies 
when compared to the dark matter halo mass function. 
The top panel of Figure [3l shows how the low-end slope 
a of the SFR function evolves as a function of red- 
shift. The evolution of dust-corrected asFR contrasts 
with the original slope of the UV LF. asFR is typically 
Aa ^ +0.13 flatter on average. Comparing aspR with 
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Figure 5. Demonstration that we can use tlie SFR function to recover tiie two inputs to the SFR function as described in Appendix A. 
We show the bolometric luminosity function at 2 ~ 4 (equivalent to the SFR function with parameters in Table [2j in the lower left panel 
and a /3 vs. bolometric luminosity relation (equivalent to a /3 vs. SFR relation) in the upper left panel. The upper right panel shows the 
/3 vs. observed UV luminosity relation we are attempting to reproduce with our simulations, and the lower right panel shows the UV LF. 
The gray area in the top left panel indicates our initial mock dataset, with each mock galaxy having a /3 drawn from a normal distribution 
with a mean < /3 > = —0.17 (M^oi + 19.5) — 1.91 and a Lbol from t he bolometri c LF in the bottom left panel. Also plotted arc the 
bolometric luminosities {red dots) and upper limits {orange arrows) of lReddv et al.l II2010I ) based on their rest-frame UV-luminositics and 
Spitzer MIPS 24^m flux measurements. The top right panel shows the simulated distribution as a function of UV luminosity, which we 
derive by dust correcting the bolometric luminosities based on the intrinsic /3 values for the individual mock galaxies. The black squares 
indicate the mean /3's we recover as a function of luminosity in these simulations, while the blue line gives the observed relation found in 
IBouwens et al.l 1)20121 ). The bottom right panel shows th at the UV LF we recover from these simulations (dot dashed lines) matches up 
well with the observed UV LF from IBouwens et al.l 1)20071 ) at 2 ~ 4 (black circles and error bars). 



similar determinations at lo wer redshift (iBothwell et al.l . 
[2011: Rcddv et al., 2008; So bral et al.l . l2012l) . we find no 
evidence for evolution with cosmic time. 

4. SUMMARY 

In this paper we combine UV LFs with dust corrections 
based on the UV-continuum slope /? and the IRX-/3 re- 
lation, to produce star formation rate functions between 
2: ~ 4 — 7. The dust correction results in a flatter faint- 
end slope and brighter M^jy. 

We find SFR functions that, at fixed (j)* , show a steady 
build-up in their value of SFR* with cosmic time, from ~ 
5 Mgyr-i at z ~ 8 to ~ 100 Mgyr-^ at 2: ~ 2, providing 
evidence that galaxies build up in a very uniform fashion 
in the first 3 Gyr (to z ~ 2). Use of the SFR function 
allows us to naturally link the evolution of normal star- 
forming galaxies at z > 4 with the population of dusty 
star- forming galaxies seen at z ~ 2 — 3. 

We thank Naveen Reddy, David Sobral and Marcel van 
Daalen for useful discussions. We are thankful to Naveen 
Reddy for sending us the information necessary to make 
a compa rison with th e brigh t z ~ 2 spectroscopic sam- 
ple from iReddv et al.l (|2010[ ). Feedback from the referee 



significantly improved this paper. We acknowledge the 
support of a NWO vrije competitie grant, ERG grant 
HIGHZ #227749, and NASA grant HST-GO-11563.01. 
PO acknowledges support from NASA through a Hubble 
Fellowship grant #51278.01 awarded by STScI. 

APPENDIX 

A. MONTE-CARLO SIMULATION TO RECOVER THE UV 
LF AND ^-LUMINOSITY RELATION 

In this study we derive SFR functions by transforming 
the UV LF using a luminosity dependent dust correction. 
However, one might suppose that the dust correction may 
be better parametrized in terms of the SFR instead of 
the UV luminosity - given that the SFR is the more 
physical quantity and dust exti nction is known to b e well- 
correlated with 'the SFR (e.g. IReddv et all . I2006f ). It is 
therefore worthwhile investigating whether we can start 
with our derived SFR functions and then recover the two 
inputs to the SFR function, i.e., the rest- frame UV LF 
and the observed /3-UV luminosity relationship. This will 
allow us to test whether the SFR functions we derive are 
sensitive to our parametrizing the mean /3 as a function 
of observed UV luminosity (instead of parametrizing it 
in terms of the SFR). We perform this test at z ~ 4, 
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where both the observed UV-luminosity function and the 
observed /3-UV luminosity relation are best determined. 

The steps in the simulation are as follows. We start 
with our z ^ 4: SFR function and we produce a mock 
data set with each galaxy in the data set having a SFR 
and (3. The initial /? distribution is drawn from a nor- 
mal distribution with a mean value equal to < l3 > ~ 
dti^ 1 (-^'^boi + 19.5) + /3j\fb„i=-i9.5 and a scatter equal to 
dint- This is shown with the grey area in the top left 
panel of Figure [5] as /3 versus bolometric luminosity and 
the simulated bolometric luminosity function is shown 
on the bottom left panel of F i gure El 

Assuming the lMeurer et al.l (jlQQgf Trelation to calculate 
dust attenuation, we obtain the observed UV luminosi- 
ties for each galaxy in this mock data set from the SFRs, 
as shown in the top right panel of Figure [S] We derive 
the P versus SFR distribution in an iterative process so 
as to best reproduce the observed /3-UV luminosity re- 
lationship with < 13 >^ -0.11 ( Muv.AB + 19.5) - 2.0 
and observed scatter CTobs — 0.34 (iBouwens et all . [20T2I ). 
The best-fit (3 vs. A/boi distribution has the form: 
< /3 > = -0.17 (A/boi + 19.5) - 1.91 and cri„t = 0.30. 

As a check on our derived (3 vs. SFR distribut i on we 
have plotted the measurements f romlReddv et al.l (I2010D 
in the top left panel of Figure El IReddv et al.l " ([20W use 
rest-frame UV and 24//m observations of z 2 galax- 
ies to study the properties of star-forming galaxies as 
a function of their bolometric luminosities. The agree- 
ment be tween the br i ght e nd of our simulated dataset 
and the IReddv et al.) (|20Tot) measurements is encourag- 
ing, though we note that at fixed bolometric luminos- 
ity one might expect z ~ 2 galaxies to have somewhat 
higher dust content and therefore redder (3 values than 
our z ~ 4 galaxies due to evolution in dust extinc t ion vs . 
SFR relation (e.g . IReddv et all . [20061 : iBuat et al.l . [2007h . 
The IReddv et al.l (|2010[ ) measurments therefore provide 
a rough upper limit to the steepness of the /3-SFR rela- 
tion at z ~ 4. 

We convert the simulated observed distribution to a lu- 
minosity function and co mpare this simulated o bserved 
UV LF with the UV LF of IBouwens eFall (pOOTi: see the 
bottom right panel of Figure [S]). The two LFs are in 
quite reasonable agreement. The characteristic luminos- 
ity M* we find is wit hin 0.13 mag of the derived LF of 
IBouwens et al.l (I2007D and t he faint-end slo p e a w e find 
is within 0.03 of the slope IBouwens et al.l (|2007f ) find. 
These differences are smaller than our current error bars 
for our SFR function parameters at z ^ 40 

^ We also considered the possibility that our results are affected 
by scatter in the IRX-/3 relation at 2 > 4. For example if we 
assume that dust extinction can be derived from the observed j3 
distribution and this j3 distribution is a simple function of the SFR 
(as discussed in Appendix A), we estimate that a 0.2 dex scatter 
in the IRX-/9 relation could possibly result in a somewhat lower 
value of SFR* (ASFR* ^ 0.09 dex). However, we could just as 
easily have assumed that dust extinction can be derived from the 
SFR of a galaxy and that scatter in the /3 distribution at a given 
SFR is due to scatter in the IRX-/3 relation which would result in 
a somewhat higher value of SFR* (which is in the opposite sense 
from the previously discussed scenario). On balance, we cannot 
correct our results for the effects of scatter in the IRX-/3 relation 
unless we have detailed knowledge about how dust extinction, /3, 
UV luminosity and SFR are correlated. In any case it would appear 
that any possible corrections are small, 0.1 dex, and one must be 
encouraged by how well our SFR, function results match up with 
the LF results a.t z 2. 



This simulation gives us confidence that our derived 
SFRs are not subject to large biases due to incomplete- 
ness at the faint end of the distribution. Furthermore, 
these simulations suggest that dust extinction is not so 
large at the bright end of the LF that we are unable to 
approximately recover the SFR function. Overall, these 
test results provide some validation for our methodology 
for performing the dust corrections. 
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